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P-1185  'A  Mathod  of  Computing  the  Inharent  Accuracy  vlth  which  a  Tlaa 
Da  lay  Can  Ba  Istlj^tad’ 

P.  Bwsrllng 

paga  6,  llna  6; 

Us lata  tha  aantanca  bag inn log  ' Tha  raaulta  for  . 'and  substitute: 

' for  tha  tlaa  being,  wa  will  also  aaauaa  that  P(t)  has  a  Pourlar 
transform  raniahing  for  |uj |  >  W.  Sararal  of  tha  statansnts  to  follow  ara 
act  rigorously  trua  if  tha  nolsa  la  band  llmitad  and  tha  function  P(t )  is 
not.  Ooa  auat  laaglna  tha  rarlous  limiting  opa rations  to  ba  carrlad  out 
in  tha  following  order:  first  tha  nolsa  bandwidth  W  is  aada  to  approaeh 
infinity  for  any  fixad  band  lialtad  function  than  tha  bandwidth  of 
adnlsslbla  funetiom  T  can  approach  infinity. ' 

Aslata  tha  two  linss  baginning  'Bacausa  of  our  assumptions . . . .  and  andii^ 
. . . .  dataminad  by1  and  substltuta: 

'Bacausa  of  our  assumptions  on  n(t)  and  F(t),  tha  sampla  spaca 
may,  in  tha  limit  as  T  — >  00  ,  ba  ragardad^  as  baing  datarminad  by  tha 
sample  spaoa  corrasponding  to' 
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SUMMARY 


Results  In  the  theory  of  statistical  estimation.  concerning  the 
greatest  lower  bound  for  the  variance  of  unbiased  estlaatora,  provide  an 
approach  to  the  problem  of  calculating  the  limits  of  accuracy  with  which 
a  time  delay  between  transmission  and  reception  of  a  waveform  can  be 
estimated. 

>  *  si^aary  is  given  of  the  requisite  results  from  estimation 
theory.  Certain  functions,  necessary  for  the  application  of  these  results 
to  the  case  of  time  delay  estimation,  are  evaluated,  assuming  the  received 
waveform  to  be  observed  against  a  background  of  additive  Gaussian  white 
noise.  A  brief  discussion  is  given  of  points  wherein  this  method  may  offer 
advantages  over  (a)  Woodward's  approach  to  the  same  problem,  and  (b )  an 
approach  based  on  the  inequality  of  Cramer-Hao. 

An  explicit  asymptotic  expression  is  calculated  for  the  minimum  error 
variance  of  unbiased  estimates  of  time  delay,  for  the  case  where  the 
a-priori  range  of  possible  tijM  delays  is  large. 
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Th a  problem  considered  Is  the  inherent  accuracy  with  which  the  ties 
delay  between  transmission  and  reception  of  a  waveform  can  be  measured. 
Woodward^ ^  considers  such  the  same  problem,  using  an  approach  based  on 
a  pesteriori  probability.  The  approach  to  be  followed  here  is  based  on 
results  in  statistical  estimation  theory  concerning  the  greatest  lower 
bound  f or  the  variance  of  estimates  of  statistical  parameters.  Before 
proceeding  to  a  more  precise  formulation  of  the  problem,  it  is  convenient 
to  susem.rlze  these  results  of  estimation  theory.  The  following  is  a 
summary,  in  a  notation  and  in  a  form  convenient  for  the  proposed  application, 
of  results  contained  in  Safe.  2,  3, and  4,  or  of  results  which  can  be 
obtained  by  straightforward  generalisation  of  these  references: 

Let  Jl  be  a  sample  space  with  points  to  ,  and  let  ^  be  a  measure 
defined  on  SX.  .  Let  TT  be  any  set  of  points  (flrlte  or  infinite)  called 
the  parameter  set,  with  individual  points  denoted  by  $  .  Let  {  p(<o,$)| 
for  $  €  IT  be  a  family  of  probability  densities  in  Jl\  with  respect  to 
the  measure  /j  .  Let  f($)  denote  a  real  valued  function  of  5  •  We  rai  i 
a  real  random  variable  ^(to)  an  unbiased  estimate  of  f ( 5 )  if 


f($) 


all  $  €  TT 


(1) 


How  pick  some  parameter  value  $„  (which  we  may  interpret  as  the  true 
value  of  |  )  and  consider 


<T 


2 

gib 


g. l.b. 


A 

n 


^(uj)  - 


f($  ) 


p(«i 


V 


(2) 
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where  g.l.b.  Beans  greatest  lower  bound  for  all  ^  satisfying  (l).  Any 

satisfying  (1)  which,  when  5  ”  SQ,  has  variance  equal  to  °2glb  [  T  \  ^  , 

is  on i led  an  unbiased  estimate  of  f($)  which  is  locally  best  for  $  «  1  • 

o 

Sow  let  us  suppose  that 


o(i  V)  = 

o 


p(<o, S  )  p(^, S' ) 


IT^TT 


dp  <  oo, 


for  all 

$0.1  ,5‘  €  TT  (3) 


(We  suppose  the  Integrand  to  be  defined  almost  everywhere  in  CL .  ) 

It  Is  easy  to  show  that  G  ($  ,  V  |  5q)  -  1  has  the  property  that 


*iaJ 


(4) 


for  any  choice  of  real  ambers  a^  and  points  S^fcTT. 

Let  us  also  denote  by  X  the  difference  between  any  two  measures  over 
TT  >  each  of  which  assigns  weight  to  only  a  finite  (but  otherwise  arbitrary) 
set  of  points  of  IT.  In  other  words,  if  f  ie  any  function  of 


jf(|>  dX($) 
TT 


i=l 


V(V 


C  5 ) 


where  a^  are  real  (positive  or  negative)  nunbers. 


Then 
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where  the  l.u.b.  seen*  lowest  upper  bound  over  el  1  possible  X  which 
assign  nan-zero  weight  to  at  least  one  point  of  TT  • 

Also,  It  can  be  shown  that  there  exists  a  sequence  such 

that  the  quantity  In  brackets  on  the  right  side  of  (6)  approaches 


Rib 


[f,  ^  ,  as  n  and  such  that 


lie 
n  —t  <*> 


G(S  1)  dXU)(  5')  -  f(S  )  -  f(  5  ) 

o  o 


TT 


Bence,  by  substitution  Into  (6), 


(7) 


gib 


|f,  50|  -  lie  j  jc(5  .5', |  s0)  dA(n)($  )  dX(n)(§') 

TT  TT 

=  He  (  Tf(|  )  -  f(5  )1  dX(n)(5  ) 

»-»«  J  L  °  J 


(8) 


It  Is  apparent  frost  (7)  that  this  sequence  |X^°^haa  the  property 

that  line  !dX^(§)  ■  0  .  (just  let  5  =  $  In  (7W- 

n-»ooj  o 

Moreover,  if  we  can  find  a  function  (or  generalized  function)  X  over  TT 
(not  necessarily  assigning  weight  to  only  a  finite  number  of  points )  such  that 
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/a(!  |  50>  dX(S')  -  f(S  )  -  f(S0> 


TT 


then,  under  certain  coal  It lone , 

2 


gib 


{f,  -  f  JG(*  *S  I  V  dX(i  }  dX(|,) 

-  J[t($  )  -  f(*o)]  dX( 5  ) 


TT 


(9) 


(10) 


Equation  (9)  is  a  generalized  Wiener-Hop?  equation.  One  auat  expect 
in  general  that  one  would  only  be  able  to  find  a  sequence  of  functions 
satisfying  the  equation  in  the  Halting  sense  of  (7),  although  there  are 
soae  cases  in  which  one  can  solve  (9)  in  closed  form.  Inability  to  solve 
(9)  exactly  should  not  trouble  one  too  ouch  in  practical  cases,  however; 
if  one  finds  any  approximate  solution  -  say  A*  -  then  a  lower  bound 
for  the  quantity 


4  satisfying  (1)  (11) 


is  obtained  by  inserting  X*  into  the  quantity  In  brackets  on  the  right 
side  of  (6).  This  is  true  because  by  (6),  the  quantity  in  brackets  on  the 
right  side  of  (6)  gives  a  lower  bound  for  the  quantity  (ll),  no  matter 
what  X  is  inserted.  Another  fact  Important  for  applications  Is  that  if 
\  is  any  (generalized)  function  satisfying  (5 ).  than  the  expression 


P-1185 

9-27-57 

5 


I  J 0<$  |  s0)  dA(S  )  d A( 5 ' ) 

TT  TT 

is  equal  to  o2glb  |  f,  for  any  f  of  the  fora 

f($  )  .  j G(5  ,!'[  SQ)  dA($')  ♦  constant. 

TT 

The  solution  A  or  \A(n)J  to  Iq.  (9)  or  (7)  will  in  general  depend 
on  SQ- 

We  will  now  proceed  to  the  precise  formulation  of  the  problem  to 
which  we  intend  to  apply  the  above  results. 

Let  F(t)  be  a  real  Valued  function  defined  over  -  1  <  00  ;  let 

a  be  a  real  number  belonging  to  same  interval  A  of  non-negative  nunbers ; 
and  let  X  be  a  real  number  belonging  to  sane  finite  interval  j_a,  b]  . 

We  suppose  the  receivec  waveform  to  be 

v(t)  «  atF(t-t)  ♦  n(t)  (12) 

where  n(t)  is  a  stationary  Gaussian  random  process  with  zero  mean  ,  and 
spectral  density  constant  and  equal  to  Nq  for  0  ^  f  1  W,  and  sero  for 

f  >W. 

We  suppose  that  v(t )  is  observed  for  -  T  < t  <  T,  where  T  Is  an  integral 
multiple  of  W  1. 

Let 

TT  *  direct  product  of  A  with  |_a,  bj 

i  -  (<*,*) 

\  ,  z  )  (  a  m  true  value  of  X  ; 

•*o  0  0  o 

T  *  true  value  o ft) 
o 


(15) 
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We  vill  investigate  the  evaluation  of  {f ,  ;  of  p rimary 

interest,  of  course,  are  the  cases 

(a)  f ( \  )  m  X  (unbiased  estimation  of  T) 

(b)  f(t)  «  «  (unbiased  estimation  of  *  ) 

Wt  ahei  1  actually  evaluate  the  results  in  the  limit  as  W  -»  oo  and 
T  —*  oo  .  (The  results  for  finite  W  and  T  are  obtained  in  the  process,  but 
become  vastly  simplified  in  the  limit  as  V  and  T  approach  infinity.  ) 
Because  of  our  assumptions  on  n(t),  the  sample  space  JT  may  be 

(«5 ) 

regarded  as  being  the  finite  dimensional  sample  space  determined  by 
e( t^)  *  v^ |  i  *  1 ,  .  .  •  ,  K ,  she  re  t  ^  ,  •••  ,  •  T  § 


and  t 
Lebesque  measure 
We  have 


Vl  “  *i 


1 

2W 


The  measure  p  vill  be  ordinary 


p(<*>,S)  -  [tVTtt]  "4*^7  fj’i  -  XT(X.±-X)  J 


(14) 


where 


wh 


(15) 


and 


0(5  1  5.) 


p(-».  S  )  p(-»»  s' ) 


p(^.  S0) 


dv1  . . .  dvK 


(16) 
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Before  proceeding  to  the  oval net Ion  of  0,  it  vould  be  worthwhile  to 
consider  am  alternative  approach  to  the  calculation  of  a  limitation  on 

inherent  accuracy  of  estimation  -  namely,  one  based  on  the  Cramer-Bao 

( 6 ) 

inequality'  For  complicity ,  we  will  consider  the  case  where  er  has 


known  value  olq,  and  t(\)  =  X  .  According  to  the  method  of  Cramer-Bao, 
a  lower  bound  for  the  quantity  (ll),  which  lower  bound  will  be  denoted  by 
0  crb  irol  ,  la  given  by  (E{ } means  expected  value  with  respect  to 


(thia  actually  holds  only  for  so-called  'regular'  estimates  ^(u>). ) 

This  lower  bound  has  the  disadvantage  of  applying  only  to  a  restricted 

class  of  estimates  ('regular');  mainly,  however,  the  disadvantage  is  that 

2  2 

there  are  some  cases  where  the  ratio  of  c  ,  to  o  is  much  lass  than 

crb  gib 

2 

unity,  so  that  a  crfa  gives  much  too crude  a  lower  bound  to  reflect 
accurately  the  true  inherent  accuracy  limitations.  We  can  illustrate 
this  by  evaluating  a2orb-  It  is  necessary  to  aasune  that  [f’(t)]2  has 
finite  integral  over  (-oo,  oo).  Inserting  (14)  into  (17),  and  passing  to 
the  limit  as  W -»  o o,  T -+&>,  gives  in  a  straightforward  manner: 


-1 


(18) 
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5cw,  there  are  many  cases  Is  which  1*  a  good  approximation  to 

2 

0  gib'  That  is  oot  •lway*  the  case  can  be  illustrated  as  follows : 

(a)  Suppose  F(t)  is  trapezoidal.  Then  a2  . — ►  0  as  the  sides  of  the 

CTD 

trapezoid  approach  the  vertical,  whereas  clearly  02  ,  would  not 

gib 

go  to  zero. 

(b  Suppose  F(t )  is  a  nodulated  carrier  of  frequency  f  .  Then 

2  ° 

0  crb^°  M  fo^*’  which  clearly  cannot  happen  for  2  „ . 

2  8*“  • 

In  this  case,  a  reflects  the  accuracy  Inherent  in  the  fine 

structure  information,  vhlch  is,  behood  a  certain  point,  over¬ 
balanced  by  ambiguity  errors. 

An  alternative  approach  to  the  question  of  inherent  accuracy  is  given 
by  Woodward^1 J  who  deals  (primarily)  with  the  case  where  F(t)  is  a  modulated 
carrier  with  bandwidth  small  compared  to  the  carrier  frequency.  He  also 
explicitly  assumes  that  fine  structure  information  is  rejected.  Using 
methods  and  criteria  baeed  on  a  posteriori  probability,  he  derives  results 


primarily  under  the  conditions 


(ii)  u(t ),  the  complex  envelope  of  F(t),  can  be  expanded  in  Taylor  series 
with  sufficiently  small  remainder  after  the  second  derivative  term. 

There  are  interesting  cases  in  which  conditions  (i)  and  (il)  do  not 
hold,  and  Woodward's  methods  seem,  from  the  computational  point  of  view, 
to  be  rather  difficult  to  extend  to  such  cases .  Also,  the  approach  based 
on  evaluating  o2^  has  two  features  which  nay  be  advantageous  in  certain 
circwstances: 

(a)  Ho  complex  functions  need  be  introduced  —  ve  deal  always  with 
the  actual  real -valued  waveform  F(t). 
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(b )  No  explicit  assiaq>ticn  need  be  — de  about  rejecting  the  fine 

■true tore  information  in  rases  where  F(t)  ia  a  Modulated  carrier. 
If  the  situation  ia  auch  that  ambiguity  error*  are  the  aain 

limitation  to  Inherent  accuracy,  thla  fact  ahould  be  autcamt i cally 

2 

reflected  in  the  value  of  a  on  the  other  hand,  if  the 

signal  energy  ia  sufficiently  great  that  fine  rtructure  error* 
are  the  aain  limitation,  thla  ahould  be  automatically  reflected. 

We  now  proceed  to  the  evaluation  of  G  by  ineertlng  {lb )  into  (l6); 
aep&rating  out  the  portion  of  the  exponent  depending  on  the  v^ ;  completing 
the  square;  and  performing  the  integration  with  respect  to  dv.^  ...  dv_. 

Hie  result  is: 


Q({  ,5'  |v  .  **P  J  P  jf  £  r2(t- 

o  i*l 


(19) 


20U  . 

X  "t<  —  W-  £  r(v  1  >  r(v  r,) 


i=l 


-2a  n 

O 


N  2W 

o  i*l 


!-  £  F<‘  -  V  F<V  T 


*2ot  o'  .  K 

~  nr  F(ti_  z0]  F(ti-  Z) 


1=1 
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W#  actually  wish  to  evaluate  the  llr.it  of  C  as  M  oo, 

gib 

T-*  a o.  W*  assume  that  this  can  be  done  by  evaluating  the  limit 
of  G  and  then  using  this  limit  in  (7)  -  (10). 

Since,  as  oo  f  At  *  — >0,  we  obtain  (assuming  F^(t)  has 

finite  integral  over  -oo,oo): 


.v|  V 


*<2 

[h(X-I')]  '  0 


[a<x-  [h(x’-  xoi] 


*/. 


(20) 


where 


2<I 


r 


oo 


F  (t)  dt 


J 

-00 


(21) 


H (X)  «  exp  |R  j>(T)J 


(22) 


with 


oo 


/>(X) 


/  r  ( t )  F(  t+X)  dt 

-00 _ 

^  00 


(23) 


F*(t)  dt 


-  ao 
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Equations  (9)  and  (10)  become 


IT 


TT  -  A  x\«,b] 

dX*(5')  -  .R  [H(f-  Tq)]  *°  dX(S') 


(24) 


and 


2 

*  gib 


dX- ( S  )  dX*(S') 


(25) 


Insertion  Into  (24)  of  f($)  =  X  deal*  with  unbiased  estimation  of 
X  i  of  f(5)  ■  oc  ,  with  unbiased  estimation  of  a  .  The  case  where  «  is 
considered  to  be  a  known  pares* ter  can  be  dealt  with  by  letting  A,  the 
a-priori  range  of  variation  of  ct  ,  be  the  single  point  0Lq.  In  this  case 
TT  -  \_a,b]  ;  Oi  ,  oc'  «  oCq  t  f($)  1  f(T)}  dX*(5)  becomes  dX*(t). 

A  class  of  lower  bounds  for  error  variance  can  be  obtained  as  follows: 
pick  any  finite  set  of  points  )  1  =  l««««n|  1°  TT  ,  and  let 
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d X •  ( % )  .  Z2  »4  J(5-  5.) 

n  1-1  1  1 


i  =  Dirac  £ -function 
in  two  dir.ensione 


real 


where  ja^  is  the  solution  to  the  set  of  simultaneous  equations 


t:  [h(t  -x ■  j] 

j-i  l  j 


otioci^2  a,/ 

°  *j  •  [H(VV]  /o  [*<»!>*  ,(V]  (27) 


Clearly  (27)  Is  a  finite  analog  of  (24). 
Aleo,  by  (6), 


{f'5o«n) 


n  r  *iVa2 

zz  [h(vtj>] 

i.3-1  J 


is  a  lower  bound  for  the  error  variance,  i.e. 


^  gib  |f’  5o;nj  “  7  gib  \  f'  *o5 


Under  certain  circuaetances,  o' 


gib  [  f ’  ^o’  n)  *  a  gib  \  * ’  M 

Cn> . 


if  the  proper  sequence  of  aets  1 5^  is  chosen.  Also,  if  any  function 

dA*($)  of  the  type  (26),  where  a.  are  any  real  numbers  whatever,  is 

1 

2  ( 

inserted  Into  (25;,  the  result  will  be  equal  to  o  jf,  for  any 

function  t(\)  of  the  fora 

r  *  *'  A2  -*  / 

f  (  $  )  *  [h(  "C  -  'c')j  0  °  dA'(S'>  ♦  constant 
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We  nay  be  able  to  chooee  dX*($)  so  that  f($)  differs  arbitrarily  little 

fran  the  desired  function.  In  sub,  even  If  one  Is  unable  to  evaluate 

exactly  the  greatest  lover  bound  for  the  variance  of,  say,  unbiased 

estimates  of  T  ,  this  aethod  enables  one  (a)  to  obtain  a  class  of  lover 

bounds  for  the  variance  which  nay  be  such  better  than,  say,  the  Craner- 

Rao  lover  bound;  or  (b )  to  obtain  the  greatest  lover  bound  for  the  variance 

of  estimates  having  a  bias  very  near  to  zero. 

It  Is  Interesting  that  (24)  can  be  solved  in  closed  fora  for  f($)  «  a  , 

and  where  X  Is  considered  known  and  equal  to  X  .  We  as suae  that  A  is 

o 

non -degenerate  and  that  Is  an  Interior  point  of  A.  Then,  (24) 
becomes 


dX*  («') 


(*  -  *  ) 
o 


(30) 


which  Is  solved  by  (since  H(o)  «=  eR) 


d  X*  (<* ) 


£•(*- 


*o> 


oC  i(eC-  «  ) 
o  o 


(31) 


where  J  and  are  respectively,  the  delta  function  and  Its  derivative; 
l.e.  for  any  function  g(«), 
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g(«  )  («-  «q)  d« 


i  (a  -  *  )  da 
o 


g(«) 


(32) 


that,  from  (25)* 

with  f($) 

«=  *  . 

—  oo  — 

n 

*  gib  if,“ol 

«o 

"  — 

N 

e 

=  T 

If2 ( t )  dt 

L-* 

75  ff2glb  ■  4- 


One  interesting  feature  of  this  result  is  that  the  answer  is 
independent  of  A,  the  a-priori  range  of  variation  of  a  ,  provided  A  is 
non -degenerate .  Iftis  Beans  that  decreasing  A  does  not  decrease  the 

error  variance  of  unbiased  estimates  of  #  —  in  other  vords,  if 
one  has  an  unbiased  estimate  attaining  the  variance  o  in  (34),  and  then 
if  A  is  decreased,  one  cannot  use  this  increase  in  a-priori  information 
to  provide  an  unbiased  estimate  of  decreased  variance.  'Siis  reflects  a 
drawback  in  this  approach  to  the  problem  of  Inherent  accuracy  —  it  does 
mot  always  adequately  reflect  the  influence  of  a-prlorl  information. 

We  may  also  obtain  in  convenient  fora  an  asymptotic  expression  for 

p 

a  as  (X  -  a)  add  (b  -  X  )  approach  infinity,  for  the  case  where  X 
gXo  o  o 

has  known  value  °C  and  f($)  «  *t  . 
o 

Tor  mathematical  convenience  we  will  assvsm  T  *  0;  later  the 

o 
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answer  for  general  T  will  bo  obtained  by  a  minor  modification  of  the 

result  for  X  -  0. 

o 

Under  the  assiswd  conditions,  we  must  sol  re  for  dX*(r)  the  equation 


H(X-t')  dV(X') 


XK(X) 


(35) 


Also, 


X  H(t)  dA*(X) 


(36) 


We  also  make  the  following  definitions: 

L(X)  .  H(t )  -  1 


of(u )  * 


.00 

-iuX 


L(X)  dX 


-  oo 


(37) 


(38) 


and  ve  assune  that 


2 

I-(X),  X  L(T) ,  X  L(X),  and  — - 

o£(u) 

are  integrable  orer  (-00,00). 


(39) 


Row  consider  the  following  dX* : 


dX’(X)  -  «(X)  dX  ♦ 


r 


jf(0) 


T  dX  -  I  dT 


(40) 
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V2  2 

b  -  a 


2^(0)  [Jf(0>  ♦  b  -  a] 


e"iuXyn(t)  dX 


i*'(u) 


*e  obtain 


H(X-T1)  dX*(x')  «  T:H(X)  ♦  remainder 


where  the  remainder  goes  to  iero  as  a  — f  -  oo ,  b  — *>  00  ,  except  for 
values  of  X  near  the  end  points  a  and  b. 

Although  this  does  not  constitute  a  rigorous  mathematical  proof, 
it  is  reasonable  to  assume  that  putting  (40)  into  (36)  will  give  the 
required  asymptotic  expression.  The  result  is 


*28i*  m  =  ?v  f  jzr  *•  ♦  ^  K  ««> «  <w*> 


-R  J  b3  -  a3 


/ . 2  2.2 

(b  -  a  ) 


3/(0)  4jC(0)  [eC(0)  ♦  b  -  a] 


(as  a  — *  -  00 ,  b  — *  00  ;  f  ( 5 )  *X  ) 


In  Yiew  of  (42),  w(x)  may  be  a  generalized  function. 


P-1185 

9-27-57 

17 


It  is  clsar  that  the  result  for  X  4  0  can  be  obtained  simply  by 

O  ’ 

substituting  b-  z  ,  a-  x  for  b,  a  respectively  in  (44).  Thus, 
o  o 


T  L(X)  dXV 


(45) 


4^(0)  [aft 0 )  -  b  -  a] 


(as  (t  -a)  — >09,  (b-T  )  — ►  oe  |  f(5)  *  X  ) 
o  o 


In  most  cases  of  interest  b  -  X  and  X  -  a  are  large  enough 

o  o 

for  the  above  approxlsations  to  hold,  while  R  is  large  enough  bo 
that  the  second  tern  in  braces  in  (45)  is  negligible. 

It  will  be  noticed  that  the  term  depending  on  [a,b]  increases 
as  (b-a)^.  On  the  other  hand,  the  estimate  a  4*  —  for  X  would  have 
bias  (X-  a-  *-  — ) ,  but  would  have  wean  square  error  not  greater  than 
(b-a)^.  Thus,  when  R  is  so  small  that  the  a-priori  range  of 
variation  of  X  is  the  main  factor  determining  mean  square  error,  the 
requirement  of  unbiasedness  is  clearly  disadvantageous.  On  the  other 
hand,  if  R  is  large  enough  so  that  mean  square  error  is  much  smaller 
than  (b-a),  as  will  be  true  in  most  cases  of  interest,  one  would 
expect  on  intuitive  grounds  that  any  optimum  estimate  based  on  a 
reasonable  criterion  would  be  approximately  unbiased  (except  when 
the  a-priori  distribution  of  X  over  [a,b]  is  known  and  non-uniform). 
The  term  in  (45)  dependent  on  [a,bj  1®  useful  chiefly  as  a 
criterion  for  how  large  R  should  be  in  order  that  the  error  variance 
be  effectively  independent  of  (b-a). 
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As  an  example,  *e  will  evaluate  <T  ^  as  given  by  (45)  for 

F(t)  «  1,  \t\i  |  t 

^  P 

»  0,  |t|>itp 


In  this  cast 


f>( X)  .  1  -  ,  |T|  ■  t 


ltl>  t 


It  la  then  readily  determined  that 

/.  \ 3  u  \  3 


^  /t  \  /t  r 

/  T2  L(T)  dr  -  4eR(/j  -  2  [r^V  <r2  ♦  2R  +  2) 


it  5 
3  p 


/(«) 


,  R  . 

2e  t 


R  £  /ut  \ 2 


2  sin  ut 


•W 


■rr-  sin  ut 

R  o 


cos  ut 


This  enables  one  to  evaluate  (45)  exactly  in  this  case. 
For  R»l,  the  expression  simplifies  to 


-R  /r/2 
•  Z  (u) 

27f  /  ^T(u) 


-  *  (»; 


(R»lj  f(S)  •=  X  5  F( t)  as  in  (46)) 
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One  final  observation  la  that  It  la  alao  poaalbla  to  daal  with  cases 
where  F(t)  Is  subjaet  to  an  unknown  dopplar  shift.  In  this  caao 

▼  (t)  «  *F£fJ(t-T)]  ♦  n(t)  (52) 

where  0  belongs  to  sobs  positive  interval  B.  In  this  case,  5  = 

The  function  G  can  be  evaluated  in  a  straightforward  wanner  slallar  to 

above,  but  cones  out  more  complicated  --It  involves  the 

function 


f(x.y) 


J F(t)  F  [x( 


—  oo 


dt 


(53) 


% 
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